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Abstract 

We study the evolution of the dark energy parameter within the scope of a spatially non-flat and isotropic 
Friedmann-Robertson- Walker (FRW) model filled with barotropic fluid and bulk viscous stresses. We have 
obtained cosmological solutions which exhibit without a big rip singularity. It is concluded that in both 
non-interacting and interacting cases non-flat open universe crosses the phantom region. We find that during 
the evolution of the universe, the equation of state (EoS) for dark energy ujd changes from ui'jj^ < — 1 to 
ui'jj^ > — 1, which is consistent with recent observations. 
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1 INTRODUCTION 

Observations of distant Supernovae (SNe la) (Perlmutter et al. 1997, 1998, 1999; Riess et al. 1998, 2000; 
Garnavich et al. 1998a, b; Schmidt et al. 1998; Tonry et al. 2003; Clocchiatti et al. 2006), fluctuation of 
cosmic microwave background radiation (CMBR) (de Bernardis et al. 1998; Hanany et al. 2000), large scale 
structure (LSS) (Spergel et al. 2003; Tegmark et al. 2004), sloan digital sky survey (SDSS) (Seljak et al. 2005; 
Adelman-McCarthy et al. 2006), Wilkinson microwave anisotropy probe (WMAP) (Bennett, et al 2003) and 
Chandra x-ray observatory (Allen et al. 2004) by means of ground and altitudinal experiments have established 
that our Universe is undergoing a late-time accelerating expansion, and we live in a priviledged spatially flat 
Universe composed of approximately 4% baryonic matter, 22% dark matter and 74% dark energy. The simplest 
candidate for dark energy is the cosmological constant. Recently, a great number of theme have been proposed 
to explain the current accelerating Universe, partly such as scalar field model, exotic equation of state (EoS), 
modified gravity, and the inhomogeneous cosmology model. There are several dark energy models which can be 
distinguished by, for instance, their EoS (w = ^) during the evolution of the universe. 

The introduction of viscosity into cosmology has been investigated from different view points (Gr0n 1990; 
Padmanabhan & Chitre 1987; Barrow 1986; Zimdahl 1996; Farzin et al. 2012). Misner (1966, 1967;) noted 
that the "measurement of the isotropy of the cosmic background radiation represents the most accurate ob- 
servational datum in cosmology" . An explanation of this isotropy was provided by showing that in large class 
of homogeneous but anisotropic universe, the anisotropy dies away rapidly. It was found that the most im- 
portant mechanism in reducing the anisotropy is neutrino viscosity at temperatures just above lO^^K (when 
the Universe was about 1 s old: cf. Zel'dovich and Novikov (Zel'dovich & Novikov 1971)). The astrophysical 
observations also indicate some evidences that cosmic media is not a perfect fluid (Jaffe et al. 2005), and the 
viscosity effect could be concerned in the evolution of the universe (Brevik & Gorbunova, 2005; Brevik et al. 
2005; Cataldo et al. 2005). On the other hand, in the standard cosmological model, if the EoS parameter 
uj is less than —1, so-called phantom, the universe shows the future finite time singularity called the Big Rip 
(Caldwell et al. 2003; Nojiri et al. 2005) or Cosmic Doomsday. Several mechanisms are proposed to prevent the 



future big rip, like by considering quantum effects terms in the action (Nojiri & Odintsov 2004; Elizalde et al. 
2004), or by including viscosity effects for the Universe evolution (Meng et al. 2007). A well known result of the 
FRW cosmological solutions, corresponding to universes filled with perfect fluid and bulk viscous stresses, is the 
possibility of violating dominant energy condition (Barrow 1987, 1988; Folomeev & Gurovich 2008; Ren & Meng 
2006; Brevikc & Gorbunovac 2005; Nojiri & Odintsov 2005). Setare (Setare 2007a, b,c) and Setare and Saridakis 
(Setare & Saridakis 2000) have studied the interacting models of dark energy in different context. Interacting new 
agegraphic viscous dark energy with varying G has been studied by Sheykhi and Setare (Sheykhi & Setare 2010). 

Recently, Amirhaslichi et al. (2011a, b); Pradhan et al. (2011); Saha et al. (2012) have studied the two-fluid 
scenario for dark energy in FRW universe in different context. Very recently Singh and Chaubey (2012) have 
studied interacting dark energy in Bianchi type I space-time. Some experimental data implied that our universe 
is not a perfectly flat universe and recent papers (Spergel et al. 2003; Bennett et al. 2003; Ichikawa et al. 2006) 
favoured a universe with spatial curvature. Setare et al. (2009) have studied the tachyon cosmology in non- 
interacting and interacting cases in non-flat FRW universe. Due to these considerations and motivations, in this 
Letter, we study the evolution of the dark energy parameter within the framework of a FRW open cosmological 
model filled with two fluids (i.e., barotropic fiuid and bulk viscous stresses). In doing so we consider both 
interacting and non-interacting cases. 

2 THE METRIC AND FIELD EQUATIONS 

We consider the spherically symmetric Friedmann-Robertson- Walker (FRW) metric as 



1 — fcr^ 



+ r\de^ Odcj)^) 



(1) 



where a{t) is the scale factor and the curvature constant k is — 1,0, +1 respectively for open, fiat and close 
models of the universe. 

The Einstein's field equations (with ^ttG — 1 and c = 1) read as 

- \r51 = -TL (2) 

where the symbols have their usual meaning and is the two-fiuid energy-momentum tensor due to bulk viscous 
dark and barotropic fluids written in the form. 

Tl ^{p + p)v^+pgl (3) 

where 

P = P- C"?, (4) 

and 

u'ui = -1, (5) 

where p is the energy density; p, the pressure; ^, the bulk- viscous coefficient; and u', the four- velocity vector of 
the distribution. Here after the semi-colon denotes covariant differentiation. 

The expansion factor 9 is defined by 6* = wfj = 3^. Hence Eq. Q leads to 

p^p- 3^H, (6) 

where H is Bubble's constant defined by 

H-'-. (7) 

a 

Now with the aid of Equations ©-((S]) and metric ([T]), the surviving field equations ^ take the explicit forms 
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and 



- + 2- + - . (9) 



Also in space-time (H]) the Bianchi identity for the bulk-viscous fluid distribution G'fj — leads to T-^ — which 
yields 

pu^ + {p + pW, (10) 



which leads to 

Using Eq. (0) in Eqs. dH) and ^ we get 
and 



p + 3i/(p + p) =0. (11) 



p=('Je-^«* + 3i/^), (12) 

p-=-(^e-2^' + 3i/^), (13) 

where p = p„i + Pd and p — p„i + po- Here p,„ and p,„ are pressure and energy density of barotropic fluid and 
P]j and Pd are pressure and energy density of dark fluid respectively. 

The equation of state (EoS) for the barotropic fluid ujm and dark field lod are given by 

Wm = , (14) 

Pm 



and 



respectively. 

From Eqs. (ITTI)-(II31) we obtain 



Now we assume 



UJD = — , (15) 

PD 



3H a? 



e-2«*. (16) 



p = a6|2 or p = 9ai^^ (17) 

where a is an arbitrary constant. Eq. ([T7|) ensure us that our universe approaches homogeneity (Collins 1977). 
This condition has also been used by Banerjee et al. (1986) for deriving a viscous-fluid cosmological model with 
Bianchi type II space time. 

Putting Eq. ([TT]) in Eq. (fT6|) and after integrating we get 

2Ht _ 3q;A^ 

which yields 

(19) 
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where A is an arbitrary constant. From Eq. ([T5|. we observe that the condition given by (fT7| restrict our study 
to the case when k = —1 (i.e. only for open universe). In the following sections we deal with two cases, (i) 
non-interacting two-fluid model and (ii) interacting two-fluid model. 



3 NON-INTERACTING TWO-FLUID MODEL 

In this section we assume that two-fluid do not interact with each other. Therefor, the general form of conser- 
vation equation (|lip leads us to write the conservation equation for the dark and barotropic fluid separately 
as, 

Pm + 3-(p„-fp™) = 0, (20) 
a 
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and 

Pd+3-{pd+Pd)^0. (21) 
a 

Integration Eq. ([20l) and using ([7|) leads to 

Prn = Poa-"^'+'^-^ or prn = poBe-'"^'+'^-^\ (22) 

where po is an integrating constant and B = ^^3(1+"™), gy using Eq. p2|) in Eqs. (|12p and (fT3|) . we first 
obtain the po and p_d in term of Hubble's constant H as 

PD - ( + 3H'^ - poi3e-3^(i+-'")*, (23) 

and 

PD = + 3i/2^ - ..™poi?e-3^(i+"-)*. (24) 

respectively. By using Eqs. ((23|) and (|24l) in Eq. (|T5|) . we can find the EoS of dark energy in term of time as 
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Figure 1: The plot of pD vs t for a — 0.1, A — 100, — 0.5 in both non-interacting and interacting two-fluid 
model 



'^^ ~ {^e-^Ht + 3^2) _ poBe-3^(i+--)* ■ 
EoS parameter for viscous DE can be written as 



(25) 



Therefore the effective 

'^^ PD " (f + 3^2) _ p„i3e-3^(i+"".)* ' ^^^^ 

The expressions for the matter-energy density flm and dark-energy density flu are given by 

4^Voi?e-i'"(iS^)(^+"-) 
''™-3if2- 31n^(3i^) ' ^''^ 

and 

„ _ po _ 6a , , iflp„Be-i"<lS''"*-> 



4 



1 1 1 1 1 1 1 1 1 f 

10 20 30 40 50 
t 

t, = Z, = 0.l ^ = 0.3 ^-0.5 



Figure 2: The plot of EoS parameter w^/ ■' vs t for po = 10, uJm — 0.5, a = 0.01, i? = 1 in non-interacting two-fluid 
model 
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Figure 3: The plot of density parameter (51) vs t ioi A = l,a = 0.01 in non-interacting two-fluid model 



respectively. Adding Eqs. (P7)) and (pS)) . we obtain 



n^n^ + no^- ^ +i. (29) 

From the right hand side of Eq. (|29)) . it is clear that for open universe, fl < 1 but at late time we see that 
— >■ 1 i.e. the flat universe scenario. This result is also compatible with the observational results. Since our 
model predicts a flat universe for large times and the present-day universe is very close to fiat, so being flat, the 
derived model is thus compatible with the observational results. 

Fig. 1 depicts the energy density of DE (po) versus t. From this figure, we observe that (pn), in both 
non-interacting and interacting cases, is a decreasing function of time and approaches a small positive value at 
late time and never go to infinity. Thus, in both cases the universe is free from big rip. 

The behavior of EoS for DE in term of cosmic time t is shown in Fig. 2. It is observed that for open uni- 
verse, the uj'j^'^ is an decreasing function of time, the rapidity of its decrease at the early stage depends on the 
larger value of bulk viscous coefficient. The EoS parameter of the DE begins in non-dark (wjj > — i) region at 
early stage and cross the phantom divide or cosmological constant {lod = — 1) region and then pass over into 
phantom {ujjj < — 1) region. The property of DE is a violation of the null energy condition (NEC) since the DE 
crosses the Phantom Divide Line (PDL), in particular depending on the direction (Rodrigues 2008; Kumar & Ya- 
dav 2011; Pradhan & Amirhashchi 2011). In theory, despite the observational constraints, extensions of general 
relativity are the prime candidate class of theories consistent with PDL crossing (Nesseris & Perivolaropoulos 
2007). On the other hand, while the current cosmological data from SN la (Supernova Legacy Survey, Gold 
Sample of Hubble Space Telescope) (Riess et al. 2004; Astier et al. 2006). CMB (WMAP, BOOMERANG) 
(Komatsu et al. 2009; MacTavish et al. 2006) and large scale structure (SDSS) (Eisenstein et al. 2005) data rule 
out that ojd ^ — 1, they mildly favour dynamically evolving DE crossing the PDL (see Rodrigues 2008; Kumar 
& Yadav 2011; Pradhan & Amirhashchi 2011; Nesseris & Perivolaropoulos 2007; Zhao et al. 2007; Coperland et 
al. 2006) for theoretical and observational status of crossing the PDL) . Thus our DE model is in good agreement 
with well established theoretical result as well as the recent observations. From Fig. 2, it is observed that in 
absence of viscosity (i.e. for ^ = 0), the universe does not cross the PDL but approaches to cosmological constant 
{ujd = —1) scenario. Thus, it clearly indicates the impact of viscosity on the evolution of the universe. 

The variation of density parameter (fi) with cosmic time t for open universe has been shown in Fig. 3. From the 
figure, it can be seen that in an open universe, is an increasing function of time and at late time, it approaches 
to the flat universe's scenario. 



4 INTERACTING TWO-FLUID MODEL 

In this section we consider the interaction between dark viscous and barotropic fluids. For this purpose we can 
write the continuity equations for barotropic and dark viscous fluids as 

Pm + 3-(pm -|-p„) = Q, (30) 
a 

and 

Pd+3-{pd+Pd)^-Q, (31) 
a 

where the quantity Q expresses the interaction between the dark components. Since we are interested in an 
energy transfer from the dark energy to dark matter, we consider Q > which ensures that the second law 
of thermodynamics is fulfllled (Pavon & Wang 2009). Here we emphasize that the continuity Eqs. ([TT]) and 
pop imply that the interaction term (Q) should be proportional to a quantity with units of inverse of time i.e 
Q oc J. Therefor, a first and natural candidate can be the Hubble factor H multiplied with the energy density. 
Following Amendola et al. (2007) and Gou et al. (2007), we consider 

Q - SHapm, (32) 

where cr is a coupling constant. Using Eq. (|32|) in Eq. ((30)) and after integrating, we obtain 

Pr,i ^ Poa ^ > or p„,^ poBe y ^ > , (33) 
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By using Eq. in Eqs. (|12l) and (IT51) . we again obtain the and in term of Hubble's constant H as 

'3/c 



Po 



_gg-3H(l+<^„-<T)t 



(34) 



and 



PD = + 3i/2^ - K, - a)poi?e-3^(i+""-'^)*, (35) 

respectively. By using Eqs. (|34|) and (|35|) in Eq. (fT5|) . we can find the EoS of dark energy in term of time as 



Again we can write the effective EoS parameter of viscous DE as 

k „-2Ht ^ 3^2^ _ 3^^ ^ ^^^^ _ ^)poSe-3H(H 
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(36) 



(37) 
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Figure 4: The plot of EoS parameter w^'^ vs t for po = 10, Wm — 0.5, a — 0.01,5 — l,a — 0.3 in interacting 
two-fluid model 

The expressions for the matter-energy density and dark-energy density VId are given by 



3ir2 



and 



D 



PD 



6a 



ln^(3^) 

respectively. Adding Eqs. ([38|) and ([39|l. we obtain 



4t2poi3e-t'"<iS^)(^+"'"-") 
31n^(3i^) 



6a 



1, 



(38) 



(39) 



(40) 



which is the same expression as in previous case of non-interacting two-fluid. Fig. 4 shows a plot of EoS 
parameter [uj'^j^^) versus t. The characteristic of lo'^j^^ in this case is the same as in the previous case. 
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5 CONCLUSION 



In this Letter, we have studied the evolution of dark energy parameter within the frame work of an open FRW 
space-time filled with barotropic and bulk viscous dark fluid. In both non-interacting and interacting cases, we 
have observed that for all values of bulk viscous coefficient, the universe has transition from non-dark region 
(u)^^ > — |) to phantom region {uj^^f < -1). In summary, we have investigated the possibility of constructing 
a two- fluid dark energy models which have the equation of state [uj'j^^) crossing - 1 by using the two-fluid 
(barotropic and bulk viscous dark fluid) naturally. Therefore, the two-fluid scenario discussed in the present 
paper is a viable candidate for dark energy. It is also worth mentioned here that in both interacting and non- 
interacting cases, our models are free from big rip. 
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